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Abstract. We discuss methods, based on the theory of vector-valued mod- 
ular forms, to determine all modular differential equations satisfied by the 
conformal characters of RCFT; these modular equations are related to the 
null vector relations of the operator algebra. Besides describing effective algo- 
rithmic procedures, we illustrate our methods on an explicit example. 



1. Introduction 

Differential equations are arguably one of the most important tools of theoret- 
ical physics. They appear in many guises, like equations of motion, conservation 
laws, etc. The usual approach to a physical theory is to deduce the governing 
differential equations starting from some basic theoretical considerations or experi- 
mental observations, and then to investigate the theory by solving these equations 
under different circumstances. But, while this is by far the most common situation, 
there are cases where one can determine the quantities of interest by some different 
method, without the knowledge of the differential equations themselves. 

Conformal characters of RCFT provide an interesting example of this phenom- 
enon. It is known [2, 22, 23] that these quantities satisfy differential equations of a 
very special kind, so-called modular equations, related to the null vector relations of 
the chiral operator algebra [7]. By analyzing the representation theory of the latter, 
one could find enough null vector relations to make the corresponding system of 
modular equations completely determined, whose unique solution should therefore 
give the conformal characters. In many examples this procedure does indeed work. 
But in some circumstances we don't have enough information about the operator 
algebra, the most extreme case being when even its existence is unknown, while 
we still have enough information to determine the would-be characters, at least 
partially, e.g. using the fact that they form a vector-valued modular form for a 
suitable automorphy factor of weight 0, which can in turn be determined from the 
would-be fusion rules. In this situation, one faces the following question: can we 
determine from the knowledge of the characters all the modular equations that they 
satisfy, and from this information infere the null vectors of the operator algebra? 
As to the second question, the precise relation of null vectors to modular equations 
has been settled in work of Gaberdiel-Keller [13] and of Zhu [28]. In this note we 
want to address the question of how to determine, for a given character vector, 
all modular differential equations satisfied by it. Our answer is based on the ma- 
chinery of vector- valued modular forms and their invariant differential operators, 
as developed in joint work with Terry Gannon [3, 4]; the relevant concepts and 
results are reviewed in the first sections. As we shall see, we can give an effective 
computational answer, which we illustrate on the well-known example of the Ising 
model. We conclude with some general comments on the applications of our results. 
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2. Vector-valued modular forms 



Recall [1, G] that the classical modular group r = SL,2(Z) of 2x2 integer matrices 
with unit determinant acts on the complex upper half-plane H = {t G C | lmr>0} 
by fractional linear transformations 



, , ar + b 

(2.1) T^7T 



cr 



d ' 



for 7= ^° *J GT. It is well known that the so-called modular curve X(T), the one- 
point compactification of the quotient H/T obtained by adjoining the cusp r = ioo, 
is a Riemann-surface of genus [18, 26]. 

In case d is a positive integer, an automorphy factor of rank d for T is a map 
q:TxH— >GLd(C) that satisfies 

(2-2) 0(7i72, t) = £»(7i,72t) e(72,r) 

for all 71 , 72 G r and tgH, and is holomorphic as a function of r; it is flat of weight 
w G K if the expression 

(2.3) e ( 7 ,r) (fe^ 1 = 0(7, r) (cr+d) - " 



is independent of r [14]. According to the above definition, a weight automorphy 
factor is nothing but a homomorphism from T to GL^C), i.e. a c?-dimensional 
matrix representation of T. From a geometric point of view, an automorphy factor 
determines a holomorphic vector bundle over the modular curve -X'(r), making 
obvious how to define direct sums and tensor products of automorphy factors. 
Note that the direct sum of two flat automorphy factors is flat only if the weights 
of the summands equal each other. 

Given an automorphy factor g of rank d, an automorphic form for g is a mero- 
morphic map X : H — > C d which satisfies the transformation rule 

(2.4) X(7T) = e( 7j r)X(T) 

for all 7GT and tgH; clearly, automorphic forms are the meromorphic sections of 
the corresponding vector bundle. An automorphic form X(r) is called weakly holo- 
morphic if it is holomorphic in the upper half-plane H, and has only a finite order 
pole at r = ioo, meaning that its Puisseux series (in terms of the local uniformizer) 
involves only finitely many negative powers; if there are no negative powers at all, 
then X(r) is holomorphic, and it is a cusp form if it vanishes in the limit r— > ioo. 
Another common appelation for weakly holomorphic forms, spread across much of 
the literature, is vector-valued modular form [4, 8, 17], emphasizing their multi- 
component nature. We shall denote by A4(g) the set of weakly holomorphic forms 
for the automorphy factor g; obviously, these form a linear space over C. 

At this point, it could be helpful to review the case of classical scalar modu- 
lar forms [1, 6, 15, 18, 20, 25]. These are automorphic forms for the rank one 
automorphy factors 

(2.5) p 2fc ( 7 ,r) = ( C T + d) 2fc 

of weight 2k, defined for any integer k. There exist holomorphic forms only for 
k>l, classical examples being the Eisenstein series [1, 25] 

Ik 

(2.6) £ 2l (,) = l--^ fft »,», 



where 

(2.7) a k (n)=J2d k 

d\ n 

is the fc-th power sum of the divisors of the integer n, and B k denotes the fc-th 
Bernoulli number. Actually, any holomorphic form may be expressed uniquely as 
a polynomial in E^t) and Eq(t). There are no cusp forms for fc<6, and there is 
a unique one (up to a multiplicative constant) for k = 6, the famous discriminant 
form 

(2-8) A(r) ^ ^ = g fi(Wf 4 • 

n— 1 

The expression of A as an infinite product shows that the discriminant form doesn't 
vanish on the upper half-plane [1, 25], and this makes it possible to construct weakly 
holomorphic forms for arbitrary k as suitable quotients of Eisensteins by powers of 
A. In particular, the so-called Hauptmodul [1, 6, 18] 



E 4 (t 



3 oo 



(2.9) J(q) = *)'' -744 = g - x +Vc(n)g" = q 1 + 196884 g + 21493760g 2 



A(r) 



=i 

is invariant under T, i.e. J{^t)~J{t) for all 7ST, is holomorphic in H, and has a 
first order pole at the cusp, hence it is a weakly holomorphic form for po ; an d every 
element of M(po) can be expressed as a univariate polynomial in J(t). Finally, we 
note that the Eisenstein series ^(r), defined by Eq.(2.6) for k = l, is equal to the 
logarithmic derivative of the discriminant form 

, . 1 d(lnA) 
(2.10) E 2 (r) = _-L-i . 

While obviously holomorphic, ^(t) is not a modular form, since it does not satisfy 
the transformation rule Eq.(2.4). 

A major simplification follows from the observation that the theory for auto- 
morphy factors of non-zero weight may be reduced to the case of zero weight via 
the so-called weight shifting trick. Indeed, thanks to the fact that the discriminant 
form doesn't vanish on the upper half-plane, it is meaningful to consider arbitrary 
fractional powers of A. This allows one to associate to the flat automorphy factor 
g of weight w the automorphy factor 

Mr) Y /l2 

of weight 0. If XGAf(g) is a weakly holomorphic form for g, then 
(2.12) X (r) = A(r)- l " /l2 X(r) 

is a weakly holomorphic form for g a , providing a one-to-one correspondence be- 
tween M.(g) and M(g ). Consequently, we can restrict our attention to forms for 
automorphy factors of weight 0, without any loss of generality. What is more, in 
many important applications the only automorphy factors that actually show up 
are of weight 0. For example, the character vector of a RCFT, formed by the genus 
one characters of the primary fields, is a vector- valued modular form[3], its weight 
automorphy factor being the celebrated modular representation determined by 
the fusion rules and conformal weights of the primary fields via Verlinde's theorem 
[27, 24]. 

For the above reasons, from now on we shall only consider flat automorphy 
factors of weight 0, i.e. matrix representations g : T — > GLd(C). We make the 
following technical assumptions on g [3, 4]: 
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(2.11) eo(7»T) = e(7,r) 



(1) £?( Q 1 j equals the identity matrix; 

(2) there exists a real diagonal matrix A, called the exponent matrix, such that 

(2.13) eQj) =cxp(27riA) . 

Note that A is not unique, since only the fractional part of its diagonal elements 
are determined by Eq.(2.13), but not the integer parts. 

If X(t) is an automorphic form for g, then (2.13) implies, taking into account 
Eq.(2.4), that the map exp(-27riAr) X(r) is periodic in r with period 1: as a con- 
sequence, it may be expanded into a Fourier series (the g-expansion of X) 

(2.14) g" A X(g) = ^X[n]g", 

where g = exp(27rir) and JL[n] 6 C d . Note that X is weakly holomorphic precisely 
when its g-expansion contains only finitely many negative powers of q. 

For any automorphy factor g, multiplication by J(t) takes the linear space Ai{g) 
to itself: in other words, AA(g) is a C[J]-module, which may be shown to be free of 
finite rank. Under suitable circumstances (that are always met in practice), there 
exists a d-by-d matrix H(r), whose columns freely generate A4(g), and such that 
the asymptotic relation 

(2.15) 3(g) q A - x as q -> 

holds for a suitable exponent matrix A. In such case we call 3(r) a fundamental 
matrix for g, and the limit 

(2.16) X = lim {q- A S(q) - q- 1 ) , 

whose existence follows from Eq.(2.15), the corresponding characteristic matrix [4]. 
As we shall see later, the numerical matrices X and A determine the fundamental 
matrix completely, and provide a useful parametrization of the different automorphy 
factors of weight 0. 

By definition, the fundamental matrix 3(r) satisfies the transformation rule 

(2.17) E(7r) = e( 7 )S(T) , 

and for each XgM(g) there exists a vector 3XsC[J] d with components that are 
polynomials in the Hauptmodul </(r), such that 

(2.18) X(t) = 3(t)0X. 

The vector 0X is called the polynomial representation of the vector- valued modular 
form X, and it will play an important role later. 

The fundamental matrix 3(r) is the most important piece of data needed to 
describe vector-valued modular forms for an automorphy factor of weight 0, so 
the question is how could one determine it. This will be achieved through the 
consideration of the invariant differential operators to be discussed in the next 
section, which are also of primary interest in the study of modular differential 
equations. 

3. Invariant differential operators 

Suppose that g is a fiat automorphy factor of rank d and weight w, and that 
XeM(g) is a vector- valued modular form for g. Except in case w = 0, the r 
derivative of X(r) fails to be a vector-valued modular form, but this can be cured 
by the introduction of a suitable correction term. Indeed, the expression 

(3.1) B w X=±f-^E 2 (r)X(r) 
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is easily shown to be a vector- valued modular form for the automorphy factor 

(3-2) 0'( 7jT ) = f? ( 7)T )fe) 

dr 

of weight w + 2. Here 

(3.3) E 2 (r) = -L^ = 1 -24,- 72^-... 

Z7TI C1T 

is the logarithmic derivative of the discriminant form A(r), cf. Eq.(2.10). 
An important consequence of (3.3) is the equivariance relation 

(3.4) T> w+ i2u (A U X) = A U D„,X , 

valid for any which allows the use of the weight shifting trick Eq.(2.12) 

discussed in the previous section; if Xe M.(g), then 

(3.5) (D.J) = D X . 

When defining higher powers of D m , one should take into account that it in- 
creases the weight by 2; as a result, one should use the recurrence relation 

(3.6) D^ +1 = B w+2n o D£ , 

Note that D™ increases the weight by 2n. To get a differential operator of order 
n that maps M.(g) to itself, one has to multiply D™ by a scalar modular form of 
weight — 2n. A suitable choice is 

(3-7) 9n(T) = A(r) n 00 , 

with nfc denoting the value of n modulo fc, and 

, n + 2n 3 + 3n 2 

(3.8) noo = . 

With the above choice of prefactors, the operators 

(3.9) V„ = c)„(t)D^ 

are invariant scalar differential operators, which means that they act on vector- 
valued modular forms component- wise, and that they map M.(g) to itself for each 
automorphy factor g. What is more, any invariant scalar differential operator may 
be expressed as a linear combination of the V n -s, with coefficients that are poly- 
nomials in the Hauptmodul 1 . In particular, this is true for the products V ra °V m . 
Some of the relevant multiplication rules read 



This follows from the observation that, for a positive integer n, any weakly holomorphic scalar 
form of weight -2n is the product of D n (T) with a form of weight 0, i.e. A4(p-2n) =Bn(i") .M(po). 
In other words, 5,i(t) generates A4(p-2n) as a C[J]-module. 

5 



VioVi 


= (J- 


-984)V 2 -± (5J+264)Vi 


V10V2 


= GM 


-744)V 3 -| (J-984)V 2 


V10V3 


= {J- 


-984)V 4 - 5 (J+744)V 3 


V10V4 


= v 5 


- i (J-984)V 4 


V10V5 


= (J~ 


-984) (J+744)V 6 - | (7J-1704)V 5 


(3.10) V 2 oVi 


= 0M 


-744)V 3 -± (5J+264)V 2 + § (J+744)Vi 


V 2 oV 2 


= GM 


-744)V 4 -§ (J+744)V 3 -§ (5J+264)V 2 


V20V3 


= v 5 - 


-(J+744)V 4 +±(J+744)V 3 




= (■/- 


-984)V 4 -± (5J+264)V 3 + § (J-984)V 2 -^ (7J-1704)V 


V30V2 


= v 5 - 


-2(J-984)V 4 + ± (5J+264)V 3 -| (J-984)V 2 


V30V3 


= (J- 


-984)V 6 -|V 5 + (J-984)V 4 - 1 3g (5J+264)V 3 . 



Of special importance are the periodicity formula 

(3.11) V„ +6 =V„oV 6 , 
which follows from D 12 A = 0, and the recursion formula 

(3.12) VioV n = a„V„+i - b„V„ , 

where a„ and b„ denote weakly holomorphic scalar modular forms of weight that 
can be expressed as the following univariate polynomials in the Hauptmodul: 



Note that, in complete accord with Eq.(3.11), the coefficients a n and b n depend 
only on the value of n modulo 6. 

Since the operators V„, together with the multiplication- by- J operator 



all map the space Ai(g) of weakly holomorphic forms to itself, Ai(g) is a module 
for the noncommutative ring IS = C[J, Vi, V 2 , ■ • ■]■ This module is necessarily of 
finite rank, since it is already of finite rank as a C[J]-module. An important result 
is that the ring D is generated by the operators J, Vi, V2 and V 3 , as a consequence 
of the relations Eqs.(3.11) and (3.10), which allow to express any operator V ra with 
n > 3 in terms of J, Vi, V2, V3, e.g. 



(3.15) V 4 = (V 2 oV 2 - V10V3 + |(J + 744)oV 3 + i(5J + 264)o V 2 ) 



It is clear from the multiplication rules Eq.(3.10) and the commutation relations 



[Vi, J]=- (J-984)o(J + 744) 
(3.16) [V 2 ,J]=- 2(J + 744)oVi + i (J + 744)o(7J-1704) 

[V 3 , J]=— 3(J- 984)o V 2 + \ (7J - 1704)o Vi - 1 (J — 984)o (7 J + 3480) , 



that the structure of the ring D does not depend on the automorphy factor: the 
ring D is universal, one and the same for all automorphy factors g. This shows that 
there is a very close connection between vector- valued modular forms (for arbitrary 
automorphy factors) and representations of D. 

Let's now turn to the determination of the fundamental matrix. Consider an 
automorphy factor g of weight 0, with fundamental matrix H(t), characteristic 



(3.13) 



0„ = (J-984)" 2 (J+74A)^^- 
b„ = '-f (J+744) + ?f (J-984) 



(3.14) 



J:M{g)^M(g) 
X(t)^J(t)X(t) 



G 



matrix X and exponent matrix A. Since each operator V„ maps A4(g) to itself, 
applying any V„ to the fundamental matrix gives a matrix whose columns are 
vector-valued modular forms: consequently, for each positive integer n there exists 
matrices T> n (r) for which 

(3.17) V„S(r) =S(t)2?„(t) , 

and whose matrix elements are weakly holomorphic scalar modular forms, hence 
polynomials in the Hauptmodul J(t); these polynomials may be determined ex- 
plicitly by comparing the g-expansions of both sides of Eq.(3.17). For example, 

(3.18) £>i(t) = (J(t)-X) (A-l) + AAf - 240 (A-l) 

(3.19) V 2 {t) = (J( r )-A , )(A-l)(A-l) + (A-i)AA , + 504(A-l)(A-§) 

and 

(3.20) 2? 3 (r) = (J(t)-X) (A— 1)(A— g)(A — |) + (A— |)(A— g)AA? 

-480(A-l)(A 2 -^f + ||) 

For any given n > 0, the relation Eq.(3.17) is nothing but a linear differential 
equation of order n that the fundamental matrix has to satisfy. A necessary and 
sufficient condition for this infinite sequence of equations to have a common solution 
is the relation 

(3.21) ViX»„(r) +X» 1 (r)X» n (r) = o n (r)X»„ +1 (r) - b„(r)D n (r) , 

which follows from Eq.(3.17) and the recursion relation Eq.(3.12); actually, only 
the n = 1 and n = 2 cases of Eq.(3.21) have to be considered, the remaining ones 
are consequences of these two. Taking into account the known expression of the 
matrices T> n (r) as polynomials in the Hauptmodul, Eq.(3.21) reduces to intricate 
algebraic relations between the matrices A and X, whose explicit form can be found 
in [4]. Provided these are satisfied, one may compute the fundamental matrix by 
solving the first order linear differential equation 

(3.22) ViS(T)=S(r)2?i(r) , 

with boundary condition Eq.(2.15): the solution will automatically solve the higher 
order equations as well. This means that S(r) is completely determined by the two 
numerical matrices A and X. 

Besides D, another ring of interest is D ho1 , spanned by linear combinations of 
powers of D with coefficients given by holomorphic scalar modular forms, i.e. suit- 
able polynomials in the Eisenstein series £4 and E§. As an algebra over C, the ring 
jjhoi j g g enera ted by the operators D, E4 and Eg, where E^ denotes multiplication 
by the Eisenstein series Ek{r). The basic commutation relations connecting these 
generators read 

[E 4 , Eb] = 

(3.23) [E 4 ,D] = ±E 6 

[E6,D] = 2E4 

Note that the elements of D ho1 usually don't preserve the weight. If an element 
of D ho1 happens to change the weight of all modular forms by 2n, then we say that 
it is homogeneous of grade n. This leads to a decomposition 

00 

(3.24) D ho1 = P^ 01 , 

n=0 
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the grade n homogeneous subspace D^ 01 consisting of the operators that can be 
written as linear combinations 

n 

(3.25) J2^-k(r)B k , 

fe=0 

whose coefficients <7fc(r) are holomorphic forms of weight 2fc; in particular, go is 
constant, while gi = 0. Each homogeneous subspace D^ 01 is finite dimensional, the 
set 

(3.26) {E°E^D"- 2a - 3b \0<a,b,n-2a-3b} 
providing a basis of it. It follows that the Hilbert series of D ho1 reads 

(3-27) (z) = (i _ z) (i _\ 2) (i _ z3) . 

4. Modular differential equations 

An interesting question, with direct relevance to physics, is to determine all 
those invariant differential operators that annihilate a given vector-valued modular 
form [12, 21, 28]. In other words, one is interested in the annihilator 

(4.1) Ann D (X) = {VeID>|VX = 0} 

of the vector-valued modular form X G M.{g)\ note that Annn(X) is a (left) ideal 
of the ring D of invariant differential operators. It is clear that, as a C[J]-module, 
Annn(X) is the inductive limit of the increasing sequence Ax CA2 C- ■ • , where A n 
denotes the set of those elements of Annn(X) whose order (as a differential operator) 

n 

does not exceed n, i.e. which can be written as a combination ^ ffc(J) with 

fe=0 

polynomial coefficients ffc£(C[J]. 

The basic idea is the following: since Xe ■M(q) implies V n X£jM(p) for all n, 
one has 

(4.2) V„X = S(r)xW 

for some Xl"! € C[ J] d , where S(r) is the fundamental matrix of g; note that X^l = JX 
is the polynomial representation of X, cf. Eq.(2.18). Let's now consider a syzygy 
S = (soi • • • ,S n ) between the vectors X^, . . . , X^, i.e. a linear relation 

n 

(4.3) 5> fc (J)X[ fe J =0 

with polynomial coefficients Sfe(J) G C[J]. Such syzygys form the syzygy module 
Syz„(X), which is finitely generated and free according to Hilbert 's syzygy theorem 
[10]. Multiplying both sides of Eq.(4.3) by H(t) from the left, and taking into 
account the definition Eq.(4.2), one gets the result that the differential operator 

n 

(4.4) V 5 = ^s fc (J)V fc 

fe=0 

annihilates X. This shows that the map s 1— > V s is a module isomorphism between 
Syz n (X) and A n . Note that this isomorphism implies that each A n is a finitely 
generated free C[J]-module. One may compute a free generating set of Syz„(X) 
using standard methods of commutative algebra, and the corresponding differential 
operators will provide a free generating set of A n . 

This is all that is needed if one is only interested in the modular differential 
equations satisfied by X up to some given order. If one is interested instead in 
the structure of the full annihilator, then one needs to understand the relation 
of A„,+i to A n , at least for large enough n. Luckily enough, this relation is fairly 
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simple. Indeed, let's denote by X n the submodule of C[ J] generated by the vectors 
X[°l , . . . , Xl"! ; clearly, these submodules form an increasing sequence Xo C Xi C • • ■ . 
But the module C[J] d is Noetherian, hence any increasing sequence of submodules 
saturates, in the sense that there exists a positive integer N (the saturation index) 
such that X„ = Xjv for all n> N, and in particular, X^l S Xn for n> N. Since 
Xat is generated by the vectors X^ , . . . , X^ , this means that for n> N there exist 
univariate polynomials p^™' GC[J] such that 

N 

(4.5) XW =^p^(J)X[ fe ] . 

k=0 

Multiplying (from the left) both sides of this equality by the fundamental matrix 
E(t), we get the equality 

N 

(4.6) V n X = ^p[ n) (J)V fc X , 

k=Q 

from which we conclude that the operator 

N 

(4.7) 2>»=V„-$>£°(.7)V fc 

k=0 

n 

belongs to A n . What is more, given an element V = fk(J) £A n , the combi- 

fc=0 

nation V — f„(J)2)„ is of order less than n, hence belongs to A n -i; it follows that 
for n> N the module A n is generated by -An-i and T) n , and the full annihilator 
Annn(X) is generated (as a C[J]-module) by the sequence iDjv+i, ■ ■ ■ and a 

generating set oi An- 

What could be said about the annihilator as a left ideal of D? Using the multi- 
plication rules Eq.(3.10), one may show that An +3 generates the full annihilator. 
This means that every modular differential equation satisfied by Xe A4(g) is a con- 
sequence of one corresponding to an element of An +3- In this respect, much more 
is true: every modular differential equation satisfied by XsA^(g) is a consequence 
of one corresponding to a generator of An+i, so that, in order to have full control 
over modular equations, it is enough to determine a generating set of the latter 
module. 

Of course, to be able to apply the above ideas, one needs first to determine 
the sequence of the X^l-s. In principle, this computation involves transcenden- 
tal operations; but, thanks to Eqs.(3.12) and (3.22), one has the simple algebraic 
recursion 

(4.8) X[" +1 I =— (p 1 (J) + b„-(J-984)(J+744)-^-lx[ ,l l . 

o„ [ aJ J 

Indeed, applying both sides of Eq.(3.12) to X, and taking into account Eqs.(4.2) 
and (3.22), one arrives at Eq.(4.8). It is straightforward to compute, starting from 
X[°l = 3X, the sequence of X^-s by using the above recursion. Once this has 
been done, the whole story boils down to some more or less elementary algebraic 
manipulations, as explained above. 

In some important applications it is not the annihilator that one is really inter- 
ested in, but rather the holomorphic annihilator 

(4.9) Ann hol (X) ={VeD ho1 I VX = 0} . 

The first thing to note is that Ann hol (X) inherits a grading from that of B) ho1 , each 
homogeneous subspace Ann hol (X) n = Ann hol (X) n B£ o1 being finite dimensional. An 
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important quantity related to this decomposition is the Hilbert-Poincare-series 



(4.10) H x (z) = J2 dim (Ann hol (X) J z n , 

n=Q 

which characterizes the rate of growth, as a function of the grade n, of the number 
of independent holomorphic operators annihilating the form X. Note that H%{z) is 
always majorized as a power series by H hoi(z) = l + z + 2z 2 + 3z 3 + - ■ ■ . 

For a given grade n > 0, an element of Ann hol (X) ra can be expressed as a lin- 
ear combination of the operators E|E^D"" 2a " 3b for < a, b, n - 2a - 3b. The 
coefficients in this linear combination satisfy a system of linear equations, whose 
coefficient matrix may be determined by considering the action of the operators 
E4EgD ,l_2a_3& on the g-expansion of X(g), and by solving this system, one gets a 
basis of Ann hol (X) n . While this direct approach is conceptually simple, its compu- 
tational complexity grows rapidly with the grade n, making it unsuitable to treat 
but the simplest cases. 

A more effective approach is based on the following observation. Any element of 
the homogeneous subspace O^ 01 of grade n can be decomposed as 

(4.11) J2 c„(«,p;e|d"- 2 »- 31 , 

a,b>0 
2a+3b<n 

where the coefficients C n (a,b) are complex numbers. If this sum annihilates the 
form X. then so does 



(4.12) C„(a, 6) D„(r)E2E*D 



n—2a — 3b 



a,b>Q 
2a+3b<r 



But the operators D n (r) E^EgD 71-20-36 don't change the weight, hence they all be- 
long to D, and being differential operators of order n — 2a — 36, they are proportional 
to V n _2a-36, the only basis element of D of that order, i.e. 

(4.13) 0„(r) E a 4 E b 6 T> n - 2a - 3b = r,„, a , 6 V„_ 2a _ 3b , 

for some prefactors 

D„(r) £ 4 (t) q E 6 ( T ) b 



(4.14) f) r 



5n-2a-3&(T) 



which are weight scalar modular forms, hence univariate polynomials in the 
Hauptmodul J(r); the precise form of these polynomials is easy to work out. Com- 
paring Eqs.(4.2) and (4.13), we get 

(4.15) K{r) E^D"- 2a " 3b X = S(r) i} n , a ,b^ ln - 2a - 3b] , 

hence the element Eq.(4.11) belongs to Ann hol (X) n if, and only if the following 
linear relation holds: 

(4.16) J2 C n (a,b)i)n,a,b^ [n - 2a - 3b] =0. 

a,b>0 
2a + 3b<n 

The sequence Xl°l,XW, . . . can be determined using the recursion relation Eq.(4.8), 
and the polynomials f) n , a ,fcGC[J] are known, so Eq.(4.16) is a linear system for the 
numerical coefficients C n (a, b): solving this system, one gets a basis of Ann hol (X) n . 
While conceptually a bit more involved, this method is much more effective than 
the direct approach based on the consideration of g-expansions. 
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5. A WORKED-OUT EXAMPLE: THE ISING MODEL 



The Ising model [7] is the Virasoro minimal model of central charge c~ 1 /2. Its 
character vector is known to be 



(5.1) X=i I f - fa 



f + fi 
f-fj 
\/2f 2 



where 

oo 

f(r)=«- 1 /«Jj( 1 + 8 »+i) , 

oo 

(5.2) fi(r) = g- 1 / 48 n(l-g" + ^) , 

n=0 

f 2 (r) = ^ 1 /24-Q (1 + gn) 

n=l 

are the classical Weber functions . These satisfy the identities 

(5.3) f! + f|=f, 

(5.4) fhh = V2 , 
and are related to the Hauptmodul trough 

(5-5) J + 744 = — A = m = m 



P 4 n 4 f 



2 



The character vector Eq.(5.1) is a modular form for the weight automorphy 
factor characterized by 



(r:) 

C:l) 



1 


1 


s/2\ 


1 


1 


i) 


V2 


-V2 




1 


-1 




1 


-1 


-V2C 3 


V2 


V2 






where £ = exp(^). A suitable exponent matrix reads 

47 

23 

2 

The fundamental matrix 



(5-6) A = 1 | 



ff + fl f 25 -fj 5 -25f-25fi 8 ( f 17 f 8_ f 24 fi _ 16f ) + i|( f 39_ f 39_ 16f 15„3 2f 15^ 

f-fl f 5 + ff-25f + 25f! 8 ( f 17 f 8 + f 24 fi _ 16f )_jL( f 39 + f 39_ 16f 15 + 32f 15) 

v ^ -(25 + fi 4 )i| f 15 f 7 (f 24_ 16) _ 16f2 4i| 
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for this automorphy factor has been determined in [3], with corresponding charac- 
teristic matrix 

/ 2325 94208 

(5.7) X = 1 275 -4096 

\ 1 -25 -23 

The character vector of the Ising model clearly equals the first column of the fun- 
damental matrix, which implies that its polynomial representation has the simple 
form 

/I 

(5.8) 3X= 

V o 

Starting from this, it is straightforward to compute the sequence of X[™l-s using the 
recursion relation Eq.(4.8), leading to the result 

, xi 2 i = ^ 



/ 23648-51J \ , / 1275J+288088 

X [3] = 856 , X [4] = -2080 

36864 V 107 ) 1769472 V -2507 

From this follows that the saturation index is N = 2, i.e. X n = X2 for n > 2, and 
that there are no syzygys between the generators of X2, i.e. Syz 2 (X) (hence A2 as 
well) is trivial. By expressing the higher X^-s in terms of X^XW and X^, one 
gets 

S3 = V 3 -^V 1 + ^(J-984) 
(5.9) 2) 4 = V 4 - ^V 2 + Jf^V! - ^(J + 744) 



Since .A3 is generated by 2)3 (because A2 is trivial), all modular equations satisfied 
by X are trivial consequences of the single equation D3X = 0. 

As to the holomorphic annihilator Ann hol (X), its Hilbert-Poincare series reads 

(5.10) i?x(z) = ^ 3 + ^ + 2z 5 + 3z 6 + • ■■= {1 _ zHl f z 2 ni _ z 3y 

showing clearly that Ann hol (X) is generated, as an ideal of D ho1 , by a single operator 
of grade 3, whose expression is (up to a multiplicative constant ) 

(5.11) 53(t)- 1 D3=D 3 -^E 4 D + ^E 6 , 

showing that, up to an irrelevant multiplicative factor, the modular equation 'S^IL — 
is a holomorphic equation. This is precisely the one that follows from the null 
vector relation for the characters of the Ising model [13]. 

6. Summary 

Differential equations satisfied by modular forms have been of interest since the 
time of Jacobi. The development of String Theory and two dimensional CFT has 
led to major advances in the application of the theory of modular forms to physics, 
and the important role of the modular equations satisfied by them has been clear 
since the early days. The work of Zhu [28] and of Gaberdiel and Keller [13] clarified 
the relation of modular equations and the structure of the operator algebra, 
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The present note addressed the question: given a vector- valued modular form X 
for some weight zero automorphy factor g, determine all modular differential equa- 
tions with (weakly) holomorphic coefficients that are satisfied by X. As we have 
seen, there exist effective algorithmic techniques, based on the general theory of 
vector-valued modular forms, that provide full control over all such modular equa- 
tions. 

An important possible use is to the existence problem of RCFT: in some ap- 
plications (e.g.[f f]), one faces the question whether RCFTs with given properties 
(fusion rules, modular properties, torus partition function, etc.) do exist or not. 
In such cases one can use the above methods to check whether there exists poten- 
tial character vectors that are consistent with the given data, have non-negative 
integral g-expansion coefficients and do satisfy suitable modular differential equa- 
tions. If no such candidate character vectors can be found, then one can conclude 
that no RCFT with a consistent operator algebra exists with the given properties; 
on the other hand, for each consistent candidate character vector the holomorphic 
modular equation that it satisfies characterize the null vector relations, hence the 
representation theory of the would-be operator algebra. 
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